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ABSTRACT 



We use Spitzer 24 //m, 70 /im and ground based Ha data for a sample of 40 
SINGS galaxies to establish a star formation rate (SFR) indicator using 70 /im 
emission for sub-galactic (~ 0.05 — 2 kpc) line-emitting regions and to investigate 
limits in application. A linear correlation between 70 /im and SFR is found and 
a star formation indicator SFR(70) is proposed for line-emitting sub-galactic 
regions as S(SFR) (M • yr" 1 • kpc~ 2 ) = 9.4 x 10~ 44 £(70) (ergs • s" 1 • kpc" 2 ), 
for regions with 12 + log(0/H) > 8.4 and S(SFR) > 10" 3 (M • yr" 1 • kpc" 2 ), 
with a 1-cr dispersion around the calibration of ~ 0.16 dex. We also discuss the 
influence of metallicity on the scatter of the data. Comparing with the SFR 
indicator at 70 /im for integrated light from galaxies, we find that there is ~ 40% 
excess 70 jum emission in galaxies, which can be attributed to stellar populations 
not involved in the current star formation activity. 

Subject headings: galaxies: ISM, ISM: structure, infrared: galaxies, infrared: ISM, 
HII regions 
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Introduction 



Interstellar dust absorbs UV and optical light and reradiates it at infrared wavelengths. 
In dusty systems, the use of UV and Ha emission to trace recent star-formation is subject to 
large uncertainties even when dust attenuation corrections are used, since these corrections 
have large scatter produced by the large range of possibl e dust content, distribution, and 



geometry relative to stars and gas present in galaxies (e .g. 



Kong et all 12004 ; 



Dale et al. 



2007 



Johnson et al. 



2007 



Meurer et a. 



Cortese et al. 



19991 



Calzetti 



2001 



20081 ). IR emission 



then becomes a reliable method to trace star formation rates (SFRs) in galaxies, where the 



UV light produced 



infrared ( IKennicutt 



py rec ent star formation is attenuated and reprocessed by dust into the 



19981 ). 



As deep galaxy surveys have often access to limited wavelength information, 
monochromatic SFR indicators offer advantages over indicators using integrated luminosity 
over extended wavelength range (e.g., FIR luminosity) by providing 'easier to use' recipes. 
Monochromatic SFR indicators based on infrared emission from both whole galaxies and 



sub-galactic regions have been investigated in deta il 



and 24 yum, by many authors (ICalzetti et al 



2005 



, particularly at the wavelength of 8 itt m 



Wu et al. 



2005 



Alonso-Herrero et al 



2006 



Perez- Gonzalez et al. 



Kennicutt et al. 



2009 



2006 



Rieke et al. 



Relano et al. 



2007 



Calzetti et al. 



2007 



Zhu et al. 



2008 



20091 ) , thanks to large samples of nearby galaxies observed 



at these wavele ngths with unprecedented resolution a nd sensitivity by the Spitzer Space 



Telescope (e.g. 



Kennicutt et al. 



2003 



Dale et al. 



20091 ) . Monochromatic SFR indicators at 



the longer Spitze r bands, 70 fim and 160 /im have also been investigated for the integrated 



light of galaxies (ICalzetti et al. 



20101 ) . These longer wavelengths may provide more reliable 
SFR indicators than either 8 um or 24 um as they are close to th e peak of dust IR emis sion 



( iRieke fc Lebofsky 



1979 



Draine et al. 



2007 



Lawton et al. 



20101 ). 



Calzetti et al. 



(120101 ) 



analyze the emission in those two bands as a SFR indicator using a sample of 189 galaxies, 
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showing that r e 



Boquien et al. 



iable SFR indicators could be established above a SFR ~ 0.1 — 0.3 M yr 1 . 



(120101 ) also presented SFR calibration of Herschel (Space Telescope) 100 /im 



and 160 jum bands for spatially resolved regions in M33. 

The launch and recent commissioning of the Herschel Space Telescope is providing 
a sensitive and high angular resolution window into the far-infrared/submm wavelength 
regime, also tracing the dust peak emission in high redshift galaxies. Upcoming ground 
facilities at mm wavelengths (ALMA and the Large Millimeter Telescope) will also provide 
increased sensitivity and coverage for deep surveys in the mm wavelength regime. Such 
surveys will probe near the rest-frame dust emission peak (~60-150 /mi), at high redshift, 
e.g., the Herschel Space Telescope will probe the rest-frame dust emission peak up to 
z~2. Thus, it is of significant interest to analyze the behavior of the Spitzer 70 /im band 
as a SFR indicator in spatially resolved sub-galactic regions, since these regions may 
resemble actively star- forming and starburst galaxies at high redshift. In this paper, we will 
investigate the 70 /im as a SFR indicator in star forming regions, with sizes from ~ 0.05 to 
2 kpc. 



A study o 



Lawton et al. 



Lawton et al. 



the 70 /iin luminosity of HII regions has been already presented by 



2010) for the Magellanic Clouds. Our analysis differs from that of 



(120101 ) in that we probe a large range of galaxies, thus matching the range 



of properties those studies that have derived SF R(70) for whole galaxies . We also use 



an un bi ased reference SFR ind icator proposed by 



Kennicutt et al 



(120071 ): 



Kennicutt et al. 



(120070 : 



Calzetti et al 



(120091 ). consisting of Ha and 24 /zm emissions. As a draw back 



our study will generally probe larger physical scales than the 



Lawton et al. 



(120 lph paper, 



although our sea 



Elmegreen et al- 



es ar e still matched to those of large star formation complexes (~100 pc, 



20061 ). 



The paper is structured as follows: §2 explains the data; §3 discusses the method for 
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deriving photometry on the sub-galactic regions; §4 compares the 70 /im against a reference 
SFR indicator and then the data are compared with a simple stellar population plus dust 
model; §5 provides the calibration of 70 fim as a SFR indicator. Discussions are presented 
in §6 and summary in §7. 

Units in this paper are, ergs • s _1 for luminosity, ergs • s _1 • kpc^ 2 for luminosity surface 
density (LSD, with symbol £), M Q • yr _1 for SFR and M Q ■ yr _1 ■ kpc -2 for star formation 
rate surface density (E(SFR) or SFRD), unless otherwise specified. We adopt a Hubble 
constant, Hq = 70 km ■ s _1 • Mpc -1 . 



Data 



2.1. Sample Selection and Description 



Our baseline sam ple is the SINGS (Spitzer Infrared Nearby Galaxies Survey, 



Kennicutt et al. 



20031 1 survey, which obtained images in both the mid-IR and far-IR with 



the Spitzer Space Telescope, plus ancillary images in the optical. For the present work, we 



are int erested in t 
in Ha. Dale et ai- 



r e sub se t of images ob t ained at 24 /im, 70 /im (with Spitzer/MIPS) and 



(120091 ): 



Calzetti et al. 



(120101 ) describe the images, including background 
subtraction for the Ha images, which were obtained from both the KPNO 2.1-m and the 
CTIO 1.5-m telescopes. Among the 75 galaxies of SINGS, we exclude all ellipticals, SO 
galaxies and some irregular galaxies, which satisfy at least one of the following conditions: 

i) only one central source could be selected but the galaxy is identified as hosting an AGN; 

ii) there are either bright star(s) across the galaxy disk, or other factors causing the quality 
of the images, when convolved to the resolution of the 70 /im images, to be degraded 

by spurious artifacts. Thus we end up with a sample of 40 galaxies; their galaxy types, 
nuclear types and adopted distances are listed in Table [TJ We don't apply any further 
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selection criterion on galaxy properties, and our 40 galaxies cover almost the full range in 
Hubble types of spiral galaxies and a few irregular galaxies, with different intensities of star 
formation, from starburst to normal star forming galaxies. Some of these galaxies have 
or may have an AGN in the center, and these central regions will be excluded from our 
analysis. What we term as 'central regions' are usually the central ~1 kpc of the galaxies, 
because of the angular resolution (~ 16") of the Spitzer 70 /im images. Hence, when 
discarding AGN-impacted nuclei, we will also be discarding regions outside of the nucleus 
that may be affected by the AGN. The large range of properties of the regions matches or 
exceed (especially in SFR surface density) the range of properties of the whole galaxies in 



our final samp 
ratios listed in 



e. The Ha imag e s are corrected for [Nil] contamination using the [Nil] /Ha 



Kennicutt et al. 



(120091 ) . 



2.2. Oxygen Abundance 



Oxygen abundances, which we will term 'me t allicit ies' in the rest of the paper, for the 



galaxies in our sample are from 



Moustakas et al. 



Moustakas et al. 



(120101 ) . They are listed in Table 8 & 9 in 



(120101 ); 1 6 out of 40 galaxi e s hay e metallicity grad ient measurements. Tw o 



calib rations were used by 



Moustakas et a. 



200J) and PT05 (IPilvugin fc Thuan 



(j2O10h . namely KK04 (IKobulnickv k Kewlev 



20051 ). to derive the metallicity information and we 



will use the average from these two methods in our analysis. Bec ause of the factor ~5 



20101), 



discrepancy in the metallicity scale resulting from the 2 methods (IMoustakas et al. 
we will avoid, to the extent possible, to refer to absolute metallicity values and mainly use 
relative values. For those with gradient information, we calculate the metallicity for each 
aperture taking into account the distance of the center of the aperture to the center of the 
galaxy, after correctin g for the projection effect using inclination information (Table 1 in 



Moustakas et al. 



20101 ). For the remaining galaxies, we assign the characteristic metallicity 
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of the whole galaxy to each aperture of this galaxy, while for 6 of the small/dwarf galaxies 



metallicity is derived from the B-band lumi nosity-met allicity 



'L-Z N 



relation and may be 



2010). Figure [2] shows the 



susceptible to additional systematics biases ( IMoustakas et al. 
distributions of the metallicities and of the metallicity uncertainties for the 597 data 
points in our sample. The uncertainties are derived from the quoted uncertainties in 



Moustakas et al. 



(120101 ) in the case of galaxies with directly measured metallicity values or 
gradients; an uncertainty of 0.2 dex is instead assigned to the metallicity value of galaxies 
derived from the L-Z relation, reflecting the factor 5 difference between the two calibration 
scales. 

We divide our sample into three sub-samples based on the adopted metallicities: a low 
metallicity sample with 12 + log(0/H) ^ 8.4 ('sub-solar', including a total of four galaxies 
and 41 apertures), an intermediate metallicity sample with 8.4 < 12 + log(0/H) < 8.8 
('solar', including a total of 25 galaxies and 425 apertures, where two galaxies' metallicities 
are from L-Z relation) and a high metallicity sample with 12 + log(0/H) 8.8 ('super-solar', 
including a total of 11 galaxies and 131 apertures, where four gal axies' metallicit i es are from 



L-Z relation). The adopted solar oxygen abundance is 8.69 from 



Asplund et al. 



J2009I). 



3. Aperture Photometry 



3.1. Source Selection 



For the comparison of mult i- wavelength images, all the 24 /xm (6") and Ha (~l"-2'' 



images have been convol ved to the sam e reso 



kernel and method from iGordon et al. 



ution as 70 //m (16") using the convolution 



( 120081 ) and registered to the same coordinate 



system and pixel size (4.5 arcsec/pix), after the global background is subtracted, which is 
determined from the mode of the pixel value distribution of the whole image. The aperture 
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size is chosen to be 16" in radius which corresponds to the FWHM of the 70 fim PSF (MIPS 
handbook 1 ) and a physical radius ~ 50 pc to 2000 pc depending on the galaxy distance 
(Tabled]). Sources are selected by manual inspection, at emission peaks of the 70 jum 
band; the other two images, at 24 /zm and in Ha, are then checked for presence of peaks in 
correspondence of the 70 /im ones. A candidate within a given aperture is accepted if it 
appears in all three images (See Fig. [Q. This will bias our sample by excluding very dust 
obscured objects (i.e. with completely absorbed Ha) and very transparent objects (with 



weak IR emission). We don't consider this a major bias in our sample, as 



Prescott et al 



(120071 ) found that in SINGS galaxies only a small fraction (~ 3%) of star forming regions is 
highly obscured. The exclusion of very transparent objects is also not considered a problem 
for our analysis, which is centered on the derivation of a SFR calibrator from IR emission, 
thus requiring presence of dust emission. Although crowding is present within our apertures, 
and often more than one HII knots are included in them, we usually can identify peaks in 
each aperture that are brighter than any other in the same aperture. We keep the overlap 
between apertures to no more than 4% of the aperture area. Some apertures in crowded 
environments need to be off-centered because of the overlap criterion and also because of 
presence of multiple emission peaks within a given aperture. The centr al regions of thes e 



galaxi es, which are classified as having or possibly having an AGN by 



Moustakas et al. 



(120 10l ). listed in Table [U are not included. With these criteria, we obtain 597 regions out of 



our sample of 40 galaxies. 

1 http:/ /ssc. spitzer.caltech.edu/mips/mipsinstrumenthandbook/MIPS_Instrument_Handbook.pdf 
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3.2. Local Background Subtraction 



Due to crowding, background annuli are difficult to determine for each 
the in fluence of a neighboring aperture. Thus, we adopt the method of 



aperture, with out 



Calzettietal 



(120051 ) to remove the local background from each aperture. Each galaxy in our sample is 
divided into several local regions usually identified as having a common environment, e. g. 
within the same spiral arm, after verifying that no sharp decrement of background, caused 
by either mosaicing problems, other data processing artifacts, or changes in the galaxy's 
environment, exists within one local background region. Then the local background for each 
aperture is determined using the mode of the pixel value distribution of the background 
region. For some more distant galaxies, only a few emission knots are resolved and are 
quite isolated. Although for these galaxies annuli around photometric apertures would 
be applicable to remove the background, we still apply our method of local background 
mode removal, for consistency with other galaxies. A comparison between our method 
and the standard background annuli method using these isolated regions shows that our 
method works within 1% accuracy in removing the background of each aperture. For 
extremely crowded regions, especially the central regions of most large spiral galaxies, 
the local background of the region is hard to determine. So higher uncertainty should be 
expected in those central regions, and extreme caution has been applied in those regions 
when determining the local background. Although the necessity of performing a local 
background subtraction is still controversial, it is not much different than a photometric 
measurement performed with an advanced annuli background to better characterize both 
the sky background and diffuse emission. We present in Appendix |A] evidence that local 
background subtraction is necessary in order to maintain a consistent behavior between low 
luminosity and high luminosity data, which is essential for performing a reliable analysis. 
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3.3. Empirically Established Aperture Correction 

We then need to define aperture corrections in order to recover the lost flux outside 
our apertures due to the significant portion of flux contained in the 70 /zm PSF wings. 
However, the aperture correction value provided by the MIPS handbook is not applicable 
to our photometry, for several reasons: i) we adopt a local background subtraction derived 
from the mode of local regions rather than annuli around apertures; ii) there is usually 
more than one emission peak within each aperture; iii) there is crowding in the aperture. 
Thus, we have to establish an aperture correction for our case. The use of a PSF at 70 /im 
to derive 'custom' aperture corrections is also not applicable, since, within each aperture, 
usually there are several emission peaks as can be seen from the high resolution images 
(see Fig. [1]), and any aperture correction will need to account for the 'extended' nature 
of our sources. Finally, each aperture suffers from contamination of the PSF wings from 
neighboring apertures because of crowding. All these make it difficult to establish the 
aperture correction simply from the theoretical PSF, as it's hard to build a reasonable 
model of PSF distribution both within and outside the aperture. We thus derive empirical 
aperture corrections using the original unconvolved high resolution Ha images, which we 
consider to represent the true flux distribution within the aperture we choose; this applies 
since the aperture size is much larger than the FWHM of the PSF of the original Ha images 
(typically l"-2"). We can then compare the original Ha images and the convolved Ha 
images to establish median aperture corrections for our aperture photometry (see Appendix 

IAD. 



3.4. Other Corrections and Error Terms 



After the aperture correction (established in Appendix |A|), the Ha pho tometry is also 
corrected for [Nil] contamination ( §2.1 j) and Galactic extinction correction (ISchlegel et al. 



negligible (IDraine 



19981 ; 10'Donnelll I1994T) . while the Galactic extinction for 24 /mi and 70 /zm is considered 



2003 



and references therein). We do not correct the Ha emissions 



for internal extinction , as we use a c ombination of Ha and 24 /zm in our analysis 



teennicutt et al. 



2007 



Calzetti et al. 



20071 ) . The Ha fluxes in this paper are the observed 



Ha fluxes as specified in these articles. 

In addition to the uncertainties from global and local backgrounds (< 2%), calibration 
uncertainties (2% for 24 /iin, 5% for 70 /zm and 10% for Ha, SINGS data release guide 2 ) and 
aperture correction uncertainties (~ 22%), we have performed other tests to determine the 
presence of possible uncertainties introduced by the single convolution kernel temperature 
used (< 1%) and misalignment or misplacing of apertures (~ 1.5%). The convolution of 
24 /im and Ha images to the 70 /im images resol ution uses a co n volut ion kernel, which 



Calzetti et al. 



( 1200(1 ) shows that there 



assumes a blackbody temperature of 50 K 3 . Since 
are two components of dust, cool (~ 20 K) and warm (~ 50 K ), we also use kernels of 25 K 
and 75 K to perform the convolution and find that less than 1% error in measured flux is 
shown between those two convolution temperature and the 50 K we use. We also shift our 
apertures by half pixel (4.5 arcsec/pixel) to see how much the photometry is changed to 
estimate the error introduced by possible misalignment or misplacing of apertures; this 
introduces at worst a ~ 1.5% difference in measured flux. 

These six error terms combined together produce our error estimate, typically ~ 24%, 
for the aperture photometry; we can easily see that the aperture correction uncertainty at 
~ 22% (Appendix E}, is the dominant source of uncertainty. 



2 http:/ /irsa.ipac. caltech.edu/data/SPITZER/SINGS/doc/sings_fifth_delivery v2.pdf 
3 http: / / dirty.as.arizona.edu/ ~kgordon/mips/conv_psfs/conv _psfs.html 
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4. Results and Analysis 

In order to establish the calibration of SFR(70), we first investigate the correlation 
between 70 /im emission and SFR and then compare the correlation with a simple model 
for dust absorption and emission of stellar light. 



4.1. Correlation Between 70 /im and SFR 



In order to determine whether the 70 /im luminosity of sub-galactic line-emitting 
regions can be used as a SFR indicator and what its limitations may be, we first need a 
reference unbiased SFR indicator for spatially resolved regions. We intentionally avoid 
TIR luminosity as the reference SFR indicator because 70 /im is a ma jor contributor to 



TIR f or most galaxies, s i nce 70 /im is near t he peak of dust e mission ( 



1979 




Draine et al. 


2007 




Lawton et al. 


2010|) 


Calzetti et al. 


( 


200 


Z); 


Kennicutt et al. 


(2009^ 



20101) . As proposed by 



Rieke fc Lebofskv 



Kennicutt et al 



fl2007h : 



(120091 ). a mixed SFR indicator, involving the 



combination of an optical and an infrared tracer of SFR, can provide an u nbiased SFR 
estima te. In this paper, we take the co mbination, L(Ha)+0. 031L(24), from 



( 120071 ) . as our reference SFR indicator. 



Calzetti et al. 



Calzetti et al. 



(120071 ) derive the coefficient 0.031 



from HI I regions. The ana 



from 



Kennicutt et al. 



ogous coefficient for the integrated light of whole galaxies, 
(j2009l ). is 0.020; the differenc e between the two is possibly due to 



2009J). Although our 



the presence of diffuse 24 /im emission in galaxies (IKennicutt et al. 
regions span nearly two orders of magnitude in size, the majority of regions are dominated 
in luminosity by HII regions or cluster s of HII regions, also on account of the fact that 



we remove the local background. The 



Calzetti et al. 



(20071) ca 



Calzetti et al. 



ibration is suitable for a 



(120101 ). and our apertures 



continuous star formation history up to 100 Myr 
are estimated to have a median crossing timescale ~ 100 Myr and at most ~ 400 Myr 
(see next section). The potential non-linearity at S(SFR) > 0.17 M Q ■ yr -1 • kpc~ 2 is an 



13 



important caveat for the application of this calibration (ICalzetti et al.ll2007l ). but the SFR of 
our a pertures reach that h igh LSD regime only in 2% of the sample. Thus the combination 



from 



Calzetti et al. 



(120071 ) is expected to be an ap propriate SFR indicato r in our case 



although we st i 



Calzetti et al. 



1 use in some cases the calibration of 



Kennicutt et al. 



(120091 ) for comparison. 



(120071 1 used Pa emission as reference SFR to calibrate the unbiased SFR 
indicator, with 33 galaxies chosen from the SINGS sample, based on the availability of 
Hubble Space Telescope (HST) Pa observations. Their sample and our sample have 19 
galaxies in common, and their sample consists mostly late type spiral and ir regular galaxies 



very s imilar to the morphology distribution of our sample. Galaxies in the 



Calzetti et al. 



( 120071 ) sample that are not included in our sample usually show centrally-concentrated 
(e.g. rings, etc.) star formation, and we had to discard them because of the lower angular 
resolution of our study coupled with the presence of central AGNs in those galaxies. One 
galaxy, NGC0024, in our s ample, is discarded by these authors due to the quality of Pa 



image ( ICalzetti et al. 



20071 ) 



The conversion used to derive SFR is then SFR(M -yr" 1 ) = 5.45 x 10 _42 L 



5.45 x 10" 42 



Calzetti et al. 



nh .(H a) + 0.031L(24)(ergs • s" 1 ) (with the stellar IMF from 



lnt-nnsi 



Kroupa 



(Ha 



2001 



20101 ). and L(X) = v L v for 24 j um an d 70 fim following the common definition 



of monochromatic flux. Adopting a 



Salpeterl (119551 ) IMF in the stellar mass range 0.1-100M o 



would increase the calibration coefficient by a factor 1.51. We use LSDs to eliminate the 
influence of the galaxy distance uncertainties, similarly to the use of LSDs in establishing 
the aperture correction. 

The relation between 70 /zm and the reference SFR indicator is shown in Figure [3j 
The linear fit in log-log space gives the correlation between 70 //m and the reference SFR 
indicator as 



log[E(70)] = (0.342 ± 0.504) + (1.036 ± 0.013)log[E(Ha) + 0.03111(24)] 



(1) 
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The error term in Equation (pO) is dominated by the dispersion of the data rather than 
the uncertainty in each data point. Although there are intrinsic differences between whole 
galaxies and sub-gala ctic regions, we a l so fit for the correlation between 70 /im and the 



Kennicutt et al. 



( 120091 ). for the integrated light of galaxies, is a 



calibration derived by 
comparison: 

log[E(70)] = (0.184 ± 0.526) + (1.042 ± 0.013)log[£(Ha) + 0.020S(24)] (2) 

This correlation is consistent with the fit in Equation ([T]) within l-a uncertainties and both 
are within ~ 3 — a of a linear correlation with a slope of unity between the two quantities. 

The slightly steeper-than-unity trend between 70 /im emission and SFR is a consequenc e 



20101 ) 



of the increased transparency of the interstellar medium at low SFR (jCalzetti et al 
At low SFR or low metallicity, the dust has lower opacity which results in a lower IR 
emission. If the fitting is constrained to high SFR, the slope asymptotically decreases to 
unity. For example, the slope is 1.004(±0.014) if fitting is constrained to log[E(70)] > 40.3 
(log[E(Har) + 0.031£(24)] > 38.9 or SFRD > 0.004 M • yr" 1 ■ kpc~ 2 ), which is consistent 
with unity. From the 70 /zm/SFR vs SFR plot (Fig. HJ), the ratio of 70 /im over SFR 
distributes almost evenly around a constant value; the l-a dispersion for the whole 
dataset is ~0.18 dex, and this value decreases to ~0.16 dex if the low metallicity galaxies 
(12 + log(0/H) ^ 8.4) are removed from the fit. In Figure [5j residuals in the 70 /zm/SFR 
ratio shows the expected correlation with metallicity (i.e. dust content) and weak or no 
correlation with L(24)/L(70) ratio (i.e. dust temperature), suggesting that dust content 
is the major contributor to the systematic scatter around the mean trend although other 
factors may still produce a small effect. We will compare our data to a simple model in 
next section to investigate the cause of the scatter in the data. 
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4.2. A Simple Model Analysis 

As can be seen from the left panel in Figure changes in metallicity produce a trend 
in the 70 /xm/SFR ratio, with higher metallicity points displaying higher 70 /zm emission in 
fixed SFR than lower metallicity ones. We can also see that the low metallicity sample data 
points, with 12 + log(0/H) < 8.4 or Z < 0.5Z Q , show a broader dispersion than the higher 
metallicity ones in Figure [3j At low metallicity, the 70 /im luminosity could be lower since 
there is not enough dust to provide the opacity to absorb UV/optical light and reradiate in 
the IR. However, other factors, such as stellar population age, dust temperature and others 
could also contribute to the dispersion. Thus, in order to further investigate the nature of 
the distribution and dispersion of t he data, we construct a simple model for dust absorption 



and emission based on the model of 



Calzetti et al. 



( 120071 ) . Details of the m odel constructio n 



are in Appendix [Bj where we use Zq — 0.0134 or 12 + log(0/H) = 8.69 (lAsplund et al. 



2009|). 



With the model, we can produce predicted 70 /im versus Ha+24 /im lines for each 
given metallicity and age. We compare the data with three models of different ages at fixed 
metallicity and vice versa, and find that the age contribution to the dispersion of data is 
almost negligible compared to the influence of metallicity. 

In Figure El three models with ages 0.01 (dotted), 0.1 and 1 Gyr (dashed), and about 
solar (~ 1.4 Z ) metallicity, are overlaid on the data points, the older the redder. If we 
consider the star formation activity within the region as triggered by a single event, the 
time required for this perturbation to propagate through out the entire region is comparable 
to the crossing time (~ 100 Myr typically for the median 700 pc physical size and from ~ 
10 Myr to 400 Myr for all the possible physical sizes, ~ 50 pc to 2 kpc, of our apertures), 
thus we consider our three age models as bracketing the likely age range of the regions in 
our sample. The difference between the three models is negligible at the low SFRD end and 
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increases towards high SFRD, but remains small at all SFRs. The small change with age is 
driven by the fact that, for constant star formation, the amount of UV and ionizing photons 
(the main contribution to the measured quantities) change negligibly for ages of 10 Myr or 
longer. Although it is not drawn on the plot, a 10 Gyr model doesn't show a significant 
difference from a 1 Gyr model but a 2 Myr model does show a significant difference from 10 
Myr model at the high SFRD end, as the ionizing population is still growing from 2 Myr 
to 10 Myr. However, we don't expect the 2 Myr model to be physically applicable in our 
sample as it is too short compared to the typical crossing time. 

In Figure we keep the age fixed at 100 Myr, of the order of a typical crossing time, 
and overlay 3 models with different metallicities, 3.7Z & , 1AZ G and 0.3Z Q , corresponding 
to the high metallicity, intermediate metallicity and low metallicity samples separately as 
the values are close to the median of these three sub-samples, on the data. The first three 
panels show the individual sub-sample overlaid with the corresponding metallicity model 
line; the models fit the average trend of data pretty well except that there is significant 
scatter in the low metallicity sample. In the last panel, these three models span a broad 
range in £(70) at fixed £(SFR), accounting at least in part for the dispersion of the data. 
This supports our argument that metallicity is the major contributor to the dispersion of 
the data. The models deviate from a linear relation; this is caused by systematic variations 
of the average stellar radiation field strength, (U, see Appendix [B]) , which increases with 
increasing luminosity and moves the dust emission peak to shorter wavelengths than 70 /im. 
The convergence of three model lines at the high luminosity end is caused by the assumption 
of a scaling between SFR and region opacity, in the sense that more active regions also 
contain more dust; thus at high luminosity our models simply probe the proportionality 
between L(70) and L(24), and no longer the scaling with metallicity. Nevertheless, each 
model is still very close to a linear relationship with unity slope. For decreasing metallicity, 
the dust infrared emission changes in several ways, which could possibly cause the large 
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scatter in the low metallicity sample. First, a decrease in metallicity will directly reduce the 
dust opacity, which reduces the extinction and thus the total amount of emission absorbed 



from stellar UV light. S econd, as the meta 



dust emission increases ( jCalzetti et al. 



2000 



l icity decreases, the effec tive temperature of 



Engelbracht et al. 



20081 ) . and thus the trend 



will start to flatten sooner. However, lacking of sufficient number of low metallicity data 
could also be an important factor of the data behavior. 

To see how the characteristics of the dust emission change as a function of SFRD, we 
plot the L(24)/L(70) ratio against SFRD, (Fig. |SJ), for both data and models. As can be 
seen from the plot, the data show a trend that is similar to that of the models, a similarity 
made even more evident when the comparison is performed with the binned data (filled 
symbols). The similarity holds at all SFR values, except for the high SFRD bin of the low 
metallicit y sample which suffe rs from low number statistics; this is in contrast with the 



results of 



Calzetti et al 



feoiOf ) who find that the L(24)/L(70) ratio flattens at hig 



h SFRD. 



Calzetti et al. 



(12010h 



We attribute the difference to the fact that the high SFRD data of 
consist mainly of LIRGs, whose dust opacity is sufficiently large that effects of self-shielding 
of the dust become important and the effective dust temperature no longer increases for 
increasing SFRD. In our case, regions at high SFRD may still display analogous properties 
of dust geometry as the low SFRD regions. 

From the comparison between the data and models, and an investigation of the 
L(24)/L(70) ratio as a function of SFR, we conclude that the observed in frared trends of 



20101 ). except for the 



the HII knots are similar to those of whole galaxies ( jCalzetti et al. 
most luminous galaxies. The differences among different HII knots appear to be mainly due 
to a luminosity scaling, with the more luminous 70 /im regions being more luminous in all 
other bands, while the systematic scatter is mostly due to differences in metallicity. This 
supports a mostly linear dependence between the 70 /im luminosity and SFR. Because the 
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low metallicity sample shows significantly larger dispersion and deviation from the mean 
trend, we will exclude the low metallicity sample to derive a SFR(70) calibration. 



5. 70 fim as A SFR Indicator 

We now derive a relation between £(70) and £(SFR) using only data with oxygen 
abundance greater than 8.4, to reduce the scatter due to the increased transparency of the 
interstellar medium. The metallicity cut excludes only 4 galaxies and 41 regions, which 
changes the relation by ~ 7.5%. For the remaining 556 regions, spanning almost 5 orders of 
magnitude in SFRD, we can approximate the trend with a unity slope relation (Fig. [9]), 

E(SFR) £(70) 



M ■ yr- 1 ■ kpc- 2 1.067(±0.017) x 10 43 erg • s" 1 ■ kpc" 2 
for 2 x 10 40 < £ 70 < 5 x 10 42 . 



(3) 



This translates into a SFR calibration, 

SFR L(70) 



(4) 



M ■ yr- 1 1.067(±0.017) x 10 43 erg • s" 1 

for 5 x 10 40 < -^(70) < 5 x 10 43 . However, the relation between luminosity and LSD has 
large dispersion due to the uncertainty in the distances. The uncertainty of the calibration 
coefficient is from the fitting and the dispersion of the data about the mean trend in 
Equations ()3]) and (@J is ~0.16 dex (dashed lines on Fig. [9]). 



Comparing our results with those of 



Calzetti et al. 



(120101 ) . we find that our calibration 



coefficient for the SFR-L(70) relation is 60% larger than theirs, which means for the 
same 70 um. luminos ity our calibration will give 60% higher SFR than the calibration of 
Calzetti et all (120101 ) . The difference can be due to the fraction of diffuse emission included 
in the measurements of whole galaxies, while we probe HH-dominated regions. Since we 
are only investigating active, star-forming regions, and we remove the local background, we 
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expect a minimal level of contamination from diffuse, non-star-forming stellar populations 
in our own analysis. 



6. Discussion 

The results and analysis in §H and $5] give a reliable mean calibration of SFR(70), 
which can be used under certain luminosity ranges and metallicity limitation (Equation [3] & 
H]). The origin and impact of the scatter around the mean correlation and the comparison 
between this calibration with other monochromatic IR SFR calibrations, are interesting 
and important issues themselves and we further discuss these issues in this section. 



6.1. Scatters in the Correlation 

From the analysis of the data and the comparison with the models, we infer that the 
systematic dispersion in the data around the mean trend is mainly due to variations in 
metallicity; variations in age and dust temperature can also produce some scatter, but 
at a much smaller level. Since metallicity and dust attenuation are correlated in first 
approximation, we should expect a relation between scatter and dust attenuation as well. 
The upper panel of Figure fTUl shows the attenuation of Ha as a function of metallicity: higher 
metallicity regions do tend indeed to have higher attenuation on average. The attenuation 
at Ha is calculated as the ratio between L nhs (Ha) and LintHnsir (Ha)=L obs (Ha)+0.031L(24), 



A(Ha) = 2.51og[L in t r i ns i C (Ha)/L obs (Ha)]((KejrmcutLel^ 



20091 ). Higher extinction results 



in more TIR emission, hence more 70 /zm emission; in fact, the lower panel of Figure fTOl 
shows that the attenuation (red squares, A(Ha) > 1; green pluses, 0.25 < A(Ha) < 1; blue 
circles, A(Ha) < 0.25) traces the scatter similar to that of metallicity. However, there may 
be objects with low extinction that have been excluded from the sample, due to our source 
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selection criterion. 



As the metallicity systematically introduces a scatter around the mean trend, we 
should expect that the calibration will change for different metallicity samples and the 
higher the sample metallicity the larger the calibration coefficient. We attempt to quantify 
this by dividing our entire sample into six different (but overlapping) sub-samples with 
different metallicity ranges and deriving the calibration for each sub-sample. The result is 
shown in Figure [TTJ The error bar on the calibration coefficient shows the dispersion in 
data of each sub-sample. From the linear fit on the figure, we could propose a metallicity 
dependent SFR calibration as 

SFR L(70) 
M Q ■ yr- 1 A(Z) x 10 43 erg ■ s" 1 U 

where A(Z) = (-8.727 ± 9.186) + (1.124 ± 1.063)(12 + log(0/H)) mean in terms of mean 
metallicity of each sub-sample. Our caveat for this calibration is the current unknown 
nature of the systematic discr epancy in the metallicity values obtained from the KK04 and 



PT05 (IMoustakas et al. 



2010) calibration scales; any change to these scales will change 



Equation ((SJ) accordingly 



r 



As a comparison, we plot on Figure [TT] the calibration constant derived by lLawton et al. 



(120101 ) for the Magellanic Clouds (with dashed line error bar). The error bar for the 



Magellanic Clouds calibration is from the dispersion of the 70 /xm/TIR r atio in their 



Lawton et al. 



f l2010h 's result, 



work. Even though our mean calibration is consistent with 
the two results are slightly inconsistent once the appropriate dependency on metallicity 
is taken into account. The discrepancy may be due to the reference SFR indicator used 



by 



Lawton et al 



(]2010l ). These authors use th e total infr a red e mission (TIR), i.e., the 



dust-absorbed starlight, and the calibration of 



Kennicutt 



(119981 ) to derive S 



ill re gions of the Magellanic Clouds. However, the SFR(TIR) as derived by 



Tis for the 



Lawton et al 



(2010) likely underestimates the true SFR in the relatively low-opacity HII regions of the 
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Magellanic Clouds as it only takes account of the obscured SFR and misses the unobscured 
SFR, thus yielding an overestimated calibration constant for SFR(70). Our test showing a 
(albeit weak) dependency of the calibration constant on the sub-sample mean metallicity 
further confirms that metallicity differences do introduce a systematic scatter in the data. 



6.2. Excess 70 /im Emission in Galaxies 



From the calibration in 



Calzettietal 



(12010( 1. a galaxy with a SFR of 1 M^yr" 1 



implies a 70 /tin luminosity of 1.725 x 10 43 ergs-s -1 . The calibration in this paper shows 
that resolved HII regions or sub-galactic star forming regions with the same total SFR of 1 
M Q -yr _1 have a total 70 /im luminosity of only 1.067 x 10 43 ergs-s -1 . The difference in these 
two calibrations reveals an average of ~40% excess 70 /im emiss ion in the galaxies. B oth 



calibrations use the Ha emission in their 'reference' SFR, and in ICalzetti et al 



(|2010| ) both 



the Ha emission and the IR emission are measured across the whole galaxy, including any 
contribution from bot h the clustered (HH) regions and the diffuse component. Thus, the 



SFR(70) calibration of 



Calzetti et al. 



(120101 ) includes contributions from both components. 
Conversely, our measurements are local, and explicitly exclude any diffuse contribution, to 
the extend possible with the angular resolution of MIPS/70; our calibration of SFR(70) 
thus also excludes any diffuse component from the galaxies. If the heating of the 70 /im 
emission in galaxies simply scaled with the Ha emission (e ither clustered or dif fuse), then 



Calzetti et al. 



:or the 



our calibration constant should be the same as that from 
same metallicity value. The presence of a significant difference between the two calibration 
constants (in the sense of 'excess' 70 /im emission in the whole galaxies) suggests that a 
portion of the 70 /tm emission from whole galaxies is in excess of what can be accounted for 
from a simple scaling of the Ha emission. Hence, we suggest that the excess 70 /tm emission 
is likely due to stellar populations that are different from those that can ionize hydrogen, 
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i.e., likely to be evolved populations older than about 10 Myr. The only other option is 
that the photons that ionize the diffuse ionized gas (DIG) can heat the dust more efficiently 
than in HII regions; we consider this scenario unlikely, as it would require a higher density 
of such photons than found in HII regions, and this is not observed. Thus, the excess 70 fim 
emission should be coming from some 'older' or 'diffuse' populations. 

Dust heated by older (>5-10 Myr), diffuse stellar populations, which are no longer 
producing ionizing photons and are not related to the most recent star formation activity, 
are still capable of heating the dust to sufficiently high temperatures that significant 
emission at 70 /im can be expected. Small star forming clusters containing only B and 
A-stars (thus, non-ionizing), but not O-stars, could also be partially responsible for the 
excess 70 /im emission observed in the integrated light of galaxies. However, we expect the 
IMF to be fully sampled when averaged over whole galaxies, and thus the effect of smaller 
star forming regions to be smoothed out. 



If we use the SFR calibration of 



Kennicutt et al. 



(120091 ) to derive a reference SFR 
(Equation [2]), we obtain a calibration coefficient in Equation (E]) of 1.285 x 10 43 , still 
suggesting that there is ~25% excess 70 /im emission in galaxies. The difference in the 
fraction also suggests a difference in L(24)/L(70) ratio between the galaxies and HII regions, 
which is discu s sed in §6.31 As our regions are only slightly larger than the ones used by 



Calzetti et al 



(120071 ) to establish the reference SFR indicator and also because of the 
application of local background subtraction, it is rea sonable to expect that we should give 
preference to the calibration of 



Calzetti et al. 



( 120071 ) for our reference SFR. In summary, 



the excess light in galaxies at 70 /im is between 25% and 40% of the total, most likely close 
to ~ 40%. 

The discrepancy between the calibration for galaxies as a whole and for resolved 
clusters of HII regions is real and significant. In Figure [T2l we compare the L(24)/L(70) 
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ratio summed up in the selected regions and that in the whole galaxy for each galaxy. 
Focusing on the AGN free galaxies (blue squares), the L(24)/L(70) ratio is systematically 
higher in the line-emitting sub-galactic sources, i.e. the active star-forming regions, than 
for the integrated light from the whole galaxy, which shows that the dust temperature in 
star forming regions is higher than that in the whole galaxy on average. A change in dust 
temperature for whole galaxies can only be driven by the presence of a stellar population 
heating the dust to a cooler temperature. We identify this stellar population as 'older' and 
'diffuse'. For the AGN contaminated galaxies (red circles) located below the one-to-one line 
in Figure H21 the dust could be heated by the central AGN to a higher temperature and 
also dominate the total IR luminosity. 

From all the analysis above, we conclude that there is at least ~25%, or more likely 
~40%, on average, excess integrated 70 /zm emission from galaxies. This excess comes from 
dust heated by older and diffuse stellar populations, which we identify as stellar populations 
not related to current star formation activity. 



6.3. Dust Temperature in HII Regions and Galaxies 



From the dust temperature versus S FR comparison (Fig. El), our data follow the models 



across the full dynamical range, while in 



Calzetti et al 



fcoiOh . a flattening on L(24)/L(70) 



occurs at the high SFR end. Their high SFR end mos tly consists of LIR Gs, whose dust 



distribution becomes optically thick even at IR bands (IRieke et al. 



20091 ) and the observed 



dust temperature starts to flatten. The difference indicates that the HII regions in our 
sample never become optically thick in the IR, even at the highest SFRDs. 
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6.4. Comparison with SFR(24) 



The 70 /im emission comprises a large portion of the TIR emission both in galaxies and 



HII regions/star-forming regions. So a legitimate 
is a better SFR indicator than other wavebands. 



question is w^ 



lether the 70 /im emission 



Lawton et al. 



(120 lOf ) shows that, in the 



Magellanic Clouds, the 70 /im emission is better than 8, 24, or 160 /im as a SFR indicator 



basec. 



on the dispersion of the data about the mean relations. Furthermore, 



Calzetti et al. 



( 120071 ) showed that the relation between the 24 /im luminosity and SFR is non-linear for 



HII regions/complexes. In the present work, we find that the 70 /im emission is linearly 
correlated with the SFR in HII regions/complexes, and shows almost a factor 2 lower 
dispersion about the mean trend than th e 24 /im emission (0.16 dex versus 0.3 dex for 



the 24 /xm emission 



Calzetti et al. 



Lawton et al. 



20071 ) . Our results thus would seem to support the 
(J2010|)'s result that the 70 /im emission is more tightly correlated with the 
SFR than the 24 /im emission in HII regions/complexes. However, we need to caution the 
reader that the reference SFR used h ere, a combination of Ha and 24 /im is different from 



the one used in 



Calzetti et al. 



(120071 ). where the extinction-corrected Pa is employed. 



For whole galaxies, iKennicutt et al.l ( 120091 ) and 



Rieke et al. 



( I2009h show that the 24 /im 



emission has a dispe rsion of only 0.12-0.16 dex about the mean trend with SFR, while 



Calzetti et al. 



( 120101 ) show that the 70 /im emission has a larger dispersion, by about 25%. 
Whether this indicates that the 24 /im emission is a better SFR indicator than the 70 /im 
emission for whole galaxies is unclear at this stage. Different calibrations rely on different 
reference SFR indicators, and there is a risk of circularity in many comparisons, both for 
whole galaxies and for HII regions/complexes. A dedicated, independent analysis using a 
consistent reference SFR indicator should be performed to solve this issue. 
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7. Summary 



A sample of 40 galaxies, with high quality Ha images and Spizter 24 /im and 70 /im 
images, has been selected from the SINGS legacy survey and 597 sub-galactic regions, 
in correspondence of peak 70 /im emission and avoiding AGN contamination, have been 
identified and measured at Ha, 24 and 70 /im. For these sub-galactic line-emitting regions 
(likely groups of HII regions), we have investigated the correlation between 70 /im dust 
emission and SFR over scales of 0.5-2 kpc to determine whether we could establish a SFR 
indicator using the monochromatic 70 /im emission. We have also investigated dependences 
on the metal content of the regions, as determined from nebular line emission, via a model 
constructed with a simple recipe for the stellar population, dust absorption, and emission. 
For the reference SFR , we have used the co mbin ation of observed 24 /im and the observed 



Ha flux calibrated in 



Calzetti et al. 



(120070 and 



Kennicutt et al. 



(120091 ). We obtain a 



relatively tight correlation between 70 /im and SFR and we provide both a mean calibration 



and a metallicity depende nt calibration. 



is similar to that found by 



Lawton et al. 



he t ight correlation between 70 /im and SFR 



( 120101 ) for HII regions in the Magellanic Clouds, 



once the difference in physical scale is taken into account. However, our detailed accounting 
of both the obscured and unobscured SFR in a large variety of galaxies and over a factor 
~10 in metallicity enables us to derive a more accurate calibr a tion c onstant between the 



Lawton et al. 



( 120 lot ). As higher resolution 



70 /im emission and the total SFR than done in 
infrared imaging will be obtained by the Herschel Space Telescope in the coming years, 
this correlation will b e further tested with higher spatial detail. Comparing with the SFR 



indicator at 70 /im of 



Calzetti et al 



(120 lOt ). which is derived for whole galaxies, we find that 
there is ~ 40% excess 70 /im emission in galaxies, which we attribute to dust heated by 
non-star-forming stellar populations. Variations in metallicity in the high and intermediate 
metallicity sample introduces a dispersion about the correlation of ~ 0.16 dex but doesn't 
affect the trend significantly. At the low metallicity end of our sample, the scatter is larger 
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as expected from the lower opacity in the regions. So for deriving the SFR calibration we 
exclude low metallicity data points. We obtain a mean calibration constant 1.067 x 10 43 
with a dispersion around the mean trend of ~ 0.16 dex. 

We avoid the regions possibly hosting or contaminated by AGN, so that the infrared 
emission in our sample is free from AGN contamination. Thus, our SFR relation will not 
be applicable to sources dominated by AGNs. Our SFR relation is established under the 
assumption of a universal stellar IMF. Adopting a different IMF will produce a different 
(scaled) calibration coefficients in Equations (j3J) & (SJ). 
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which is operated by the Jet Propulsion Laboratory, California Institute of Technology 
under a contract with NASA. This work has been partially supported by the JPL, Caltech, 
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Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and Space Administration. 
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A. Empirical Aperture Correction 

As stated in §3.3[ we use photometry on unconvolved high resolution Ha images to 
establish our empirical mean aperture corrections for the photometry in 16" apertures on 
the 24 /im and 70 /xm images. 
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Figure [13] shows how the difference between the original unconvolved photometry and 
convolved photometry is distributed as a function of the original photometry of Ha as a 
reference (with local background subtraction). If we had one perfect point source in each 
of our apertures, the expected correction to the flux would be 0.3 dex (from the MIPS 
handbook, using background annuli from 18" to 39"). Most of the data distribute a little 
below the horizontal line of ~ 0.3 dex; this is expected since a method for local background 
removal ( §3.2[) makes the aperture correction smaller as it receives a smaller contribution 
from the PSF wings than the background annuli method used in the MIPS handbook. For 
decreasing luminosity we can also see that the difference tends to be smaller or even below 0, 
which means the convolved photometry is getting larger than the original one. This trend is 
due to the contamination of neighboring apertures, and the fainter the aperture the higher 
the contamination. For the extremely faint apertures, the trend flattens, as these regions 
can only be 'identified' in relatively uncrowded regions. Since within a given aperture there 
is more than one emission peak for most cases, most apertures are not perfectly centered on 
the emission peak. This will produce a dispersion in the aperture correction due to different 
distribution of emission peaks within a given aperture. Those few data points indicating 
aperture corrections greater than a factor ~ 2.5 (Fig. [TBI are due to the relatively faint 
emission within the regions, because either the S/N is low or the aperture is substantially 
off-centered to avoid overlapping with adjacent bright sources. This is also reflected in the 
large error bars in the photometry of these data points. In summary, the trend of Figure [13] 
indicates that in addition to the overall constant aperture correction due to the PSF wings 
loss outside our selected apertures (observed in the high luminosity region of the plot) there 
is also a surface brightness dependent aperture correction, indicative of the contamination 
effects of neighboring apertures. 

Figure dH shows the same plot without local background subtraction for either 
photometric measurement. The significant difference in this figure relative to Figure [13] 
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is the flaring of the data points at the low end of the luminosity distribution, showing a 
markedly different behavior from the high end. Furthermore, the low-luminosity end is 
still more than 0.1 dex, on average, below the mean of the high- luminosity end in Figure 
mi again showing evidence for contamination from neighboring sources. We interpret 
the differences and similarities between Figure [13] and [14] as indicating the necessity of 
subtracting a local background from all apertures and the photometry measurements used 
hereafter are all with local background subtraction. 

We finally test our aperture corrections for an additional source of bias: the impact 
of sources compactness. We can technically argue that a more compact source will show a 
larger aperture correction than a more diffuse source, when convolved to lower resolutions. 
Thus, our 'decrease' in aperture correction with decreasing luminosity may simply indicate 
that less luminous sources are more diffuse than more luminous sources. To test for this 
effect we identify within our apertures all those that are dominated by a single source, i.e. 
the sources represent more than half of the flux in the entire aperture, and measure their 
FWHM in the original unconvolved images (58 sources). We also identify an additional 96 
apertures where two sources dominated the flux but one is significantly brighter than the 
other. Figure [15] shows the first group of 58 sources as filled circles, and the second group 
as empty circles. We can see that both groups do not show a dependence of the aperture 
correction on FWHM or the compactness of the sources within each aperture. A linear fit 
(dashed line) through the 58 filled circles, the clean sample, shows consistency with a slope 
equal to zero. We thus conclude that our original interpretation, the decreasing aperture 
correction with decreasing luminosity indicates increasing contamination from neighboring 
sources, is the most likely. We will then use luminosity surface density to quantify aperture 
corrections for our data. 

Based on the discussions above, we model our aperture correction as two parts, a 
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constant aperture correction for recovering the loss of PSF wings and a surface brightness 
dependent aperture correction for removing the contamination of PSF wings from 
neighboring apertures. 

A.l. Constant Aperture Correction 

Since the most luminous apertures are proportionally less affected by contamination 
from neighboring apertures, we use the high end of Figure [T3] as a gauge to establish the 
constant aperture correction. From Figure [13] we can see that the 17 high LSD apertures, 
log[S orig (Ha)] > 40, more or less lie along one horizontal line with some dispersion. So we 
take the low end, 39.971, of these 17 points, ~ 3% of the total data (blue points on Figure 
[TB"]) . as the cut off LSD and we take the mean aperture correction of these points and get 
an aperture correction as 1.788. The choice of LSD cut-off is arbitrary, but it won't change 
the final correction much as long as it stays in the high luminosity end. 

A. 2. Surface Brightness Dependent Aperture Correction 

After the constant aperture correction is applied, Figure [16] shows the difference 
between original and c-corr (constant-corrected, see section lA.lj) photometry versus c-corr 
photometry. We change the x-axis from original photometry to the c-corr photometry 
because when we apply the surface brightness dependent aperture correction on the 24 /im 
and 70 /im data, we need to rely on the c-corr photometry as we have no access to the 
original high resolution photometry. We fit the data with log[S orig (Ha)] < 39.971 with a 
linear function giving a best fit: 

log[S orig (Ha)] - log[S c _ corr (Ha)] = 0.103[log[S c _ corr (Ha)] - 39.971] 



-30 - 



as indicated by the black line in Figure [161 Thus we have determined the empirical 
aperture correction for Ha as two steps, 

1) the constant aperture correction 

£ c _ corr (#a) = 1.788 x ^ conv (Ha) (Al) 
to correct the convolved photometry to the c-corr photometry and 

2) the surface brightness dependent aperture correction 

^ corr {Ha) = io - 103 ^— (Ha)]-39.97i] x Y, c ^ corr {Ha) (A2) 

for log[£ c _ COIT (Ha)] < 39.971 to correct the c-corr photometry to 1-corr (luminosity- 
corrected) photometry. 

Figure [T7] shows the corrected (1-corr) photometry and the difference between the 
corrected (1-corr) and original photometry versus the original photometry. We estimate 
the I — a aperture correction uncertainty to be ~ 22.3% by calculating the RMS of the 
deviation from the unity slope line (lower panel on Fig. IT71) . 

We further test the robustness of an aperture correction by considering only the 
brightest half of our sample (in terms of LSD) and repeating the above analysis. The 
results are shown in Figure [181 with the brightest half of the sample shown as circles. As 
can be seen from the Figure, the final results do not change whether the entire sample or 
the brightest half is used to establish the aperture correction. 

We then apply these aperture corrections to 70 fim and 24 /jm, according to a relation 



that 0.031 x £(24) ~ L(Ha) and 0.2 x L(70) ~ L(24) f lCalzetti et al.l 120071 . 120101 ): an 



analysis of the distribution of 70 fim and 24 fim versus Ha photometry on our data gives 
similar factors. The constant aperture correction remains the same for these two IR bands 
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while the second step, Equation flA2j) . becomes 

S ; _ corr (24) = l 0- 1 03[^[0-03 12 c- C o rr (24)F39.971] x Y, c _ corr {2A) (A3) 

for log[0.031£ c _ corr (24)] < 39.971, and 

S ,_ corT (70) = l 0- 1 03[log[0.2x0.031 Sc _ co „(70)]-39.971] x ^^(70) ( A 4) 

for log[0.031 x 0.2£ c _ corr (70)] < 39.971. 

As we have no information about the original photometry of 70 /iin and 24 /im, the 
1-corr photometry will be a proxy to represent the true photometry of these two bands, 
and for consistency, we will also use the 1-corr photometry of Ha in the following analysis. 
All the luminosity and LSD values hereafter are derived from the 1-corr photometry unless 
otherwise specified. 



B. Construction of the Simple Model 



The simple model is constructed based on the model of lCalzetti et al.l (120071). and it h as 



three basic ingred ients: stellar popula tion models (STARBURST99, 



Leitherer et al. 



dust attenuation (Calzetti et al. 



(SED) (models from 



2000) 



Draine fc Li 



19991 ). 



2003). 



and dust emission spectral energy distribution 



For the stellar population models, we adopt a continuous star formation history based 
on the analysis of typical crossing time scale for our regions, calculated by using the sound 
speed and the aperture size, which is ~ 100 Myr. We take this time scale as representative 
of our stellar population, but, for completeness, we generate constant star formation 
models between 2 Myr and 10 Gyr. For the stellar atmosphere metallicity provided by 
the STARBURST99 code, we choose Z=0.004, 0.02 and 0.05, whic h is roughly 3Z^ 



1.4Z and 3.7Z Q , where Z Q is equivalent to 12 + log(0/H) = 8.69 (lAsplund et al 



2009J) 
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We derive Ha luminosities by integrating over the ionizing photon of the spectrum, and 
assuming case B recombination. After this step, we have a stellar population SED and the 
intrinsic Ha luminosity. 



Then we take the population SED and apply the attenuation curve by ICalzetti et al. 



( 20001 ) to get the TIR emission, assuming all the absorbed UV-optical light has been 



reradiated in the IR. The E(B-V) value needed for the at tenuation cury e is ca lculated from 



the ionizing photon number using an empirical relatio n in 



Calzetti et al. 



of E(B-V) with metallicity is also taken into account (ICalzetti et al. 



(120071 ) . The scaling 



20071 ). For the nebular 



line attenuation s, we adopt dif f erenti al extinction with the stellar continuum (E S =0.44E S 



as presented in 



Calzetti et al. 



(120001 ) plus the Milky Way extinction curve. With the 



attenuation curve applied, we have the TIR luminosity and the observed Ha luminosity. 



We finally use the dust model from lDraine fc Lil ( 120071 ) to determine the L24/TIR and 
L70/TIR fractions and to further get 24 fim luminosity and 70 /im luminosity. The dust 
emission models are parametrized as a function of qpAH, the fraction of PAH molecules, 
and U, the average stellar radiation field strength. Qpah has only a modest impact on 
our output, and we adopt two values: qpAH = 4.6% for solar and super-solar metallicity 
models, and qpAH = 0.47% for the sub-solar metallicity model. We obtain estimates of U 
by integrating the st ellar population SED using an approach similar to that presented in 



Calzetti et al 



(120071 ) . By combining these inputs with the models, we have the observed 



24 /im luminosity and 70 /im luminosity. 

Finally, we assume each model is observed within an unresolved 700pc aperture, which 
is the median physical size of all the apertures, and get the LSD values for Ha, 24 fim and 
70 fim. of the model, to compare with the data. The change in the adopted aperture for 
deriving model LSDs only causes a small shift (smaller than the dispersion of data) of the 
model lines. Also a test on the data of two fixed physical sizes (~ 250 pc and ~ 650 pc 
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with ~ 100 apertures each) reveals that the correlations for data of difference physical sizes 
are consistent with each other. 
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Table 1. Sample 



Galaxy Type Nuc D Size Oxygen Abundance # Regions Sub-sample 

(Mpc) (pc/16") (Information) 
(1) (2) (3) (4) (5) (6) (7) (8) 



NGC0024 


SAc 


SF 


8.2 


636 


1 


3 


M 


NGC0337 


SBd 


SF 


24.7 


1916 


1 


5 


M 


NGC0628 


SAc 


SF 


11.4 


884 


2 


40 


M 


NGC0925 


SABd 


AGN 


10.4 


807 


2 


23 


M 


NGC1097 


SBb 


AGN 


16.9 


1311 


2 


20 


M 


NGC1512 


SBab 


AGN 


10.4 


807 


1 


7 


H 


NGC1566 


SABbc 


AGN 


18 


1396 


3 


15 


H 


NGC1705 


Am 


SF 


5.8 


450 


1 


2 


L 


NGC2403 


SABcd 


SF 


3.5 


271 


2 


41 


M 


Ho II 


Im 


SF 


3.5 


271 


1 


8 


L 


NGC2798 


SBa 


AGN 


24.7 


1916 


1 


1 


M 


NGC2841 


SAb 


AGN 


9.8 


760 


2 


9 


H 


NGC2976 


SAc 


SF 


3.5 


271 


1 


9 


M 


NGC3049 


SBab 


SF 


19.6 


1526 


1 


1 


H 


NGC3190 


SAap 


AGN 


17.4 


1350 


3 


3 


H 


NGC3184 


SABcd 


SF 


8.6 


667 


2 


27 


H 


IC2574 


SABm 


SF 


3.5 


271 


1 


10 


L 


Mrk33 


Im 


SF 


21.7 


1683 


1 


1 


M 


NGC3351 


SBb 


SF 


9.3 


721 


2 


8 


H 
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Table 1 — Continued 



Galaxy Type Nuc D Size Oxygen Abundance # Regions Sub-sample 

(Mpc) (pc/16") (Information) 
(1) (2) (3) (4) (5) (6) (7) (8) 



NGC3521 


SABbc 


AGN 


9 


698 


2 


22 


M 


NGC3627 


SABb 


AGN 


8.9 


690 


1 


8 


M 


NGC3938 


SAc 


SF 


12.2 


946 


3 


22 


M 


NGC4254 


SAc 


SF 


20 


1551 


2 


20 


M 


NGC4321 


SABbc 


AGN 


20 


1551 


2 


24 


H 


NGC4450 


SAab 


AGN 


20 


1551 


3 


5 


H 


NGC4536 


SABbc 


SF/AGN 


25 


1939 


1 


8 


M 


NGC4559 


SABcd 


SF 


11.6 


900 


2 


17 


M 


NGC4579 


SABb 


AGN 


20 


1551 


3 


7 


H 


NGC4631 


SBd 


SF 


9 


698 


1 


14 


M 


NGC4725 


SABab 


AGN 


17.1 


1326 


1 


22 


M 


NGC5055 


SAbc 


AGN 


8.2 


636 


2 


24 


M 


NGC5194 


SABbc 


AGN 


8.2 


636 


2 


35 


H 


NGC5474 


SAcd 


SF/AGN 


6.9 


535 


1 


6 


M 


NGC5713 


SABbcp 


SF 


26.6 


2063 


1 


1 


M 


IC4710 


SBm 


SF 


8.5 


659 


3 


8 


M 


NGC6822 


IBm 


SF 


0.6 


47 


1 


18 


L 


NGC6946 


SABcd 


SF 


5.5 


427 


2 


37 


M 


NGC7331 


SAb 


AGN 


15.7 


1218 


2 


22 


M 
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Table 1 — Continued 



Galaxy Type Nuc D Size Oxygen Abundance # Regions Sub-sample 

(Mpc) (pc/16") (Information) 
(1) (2) (3) (4) (5) (6) (7) (8) 



NGC7552 SAc SF 22.3 1730 
NGC7793 SAd SF 3.2 248 



1 

2 



3 

41 



M 
M 



References. — (1) Kennicutt et all (120031 ): (2) Moustakas et all (J2010j). 



Note. - Col. (1) SINGS galaxy nam e; Col. (2) Morpholog ical type; Col. (3) Adopted 



nuclear optical s pectral classificati o n from 



Moustakas et al. 



it is SF/AGN in 
are adopted from 



Moustakas et al. 



Kennicutt et al 



(12010 



( 120101 ): type AGN is adopted when 
; Col. (4) Distance : morphological type an d distance 



(120031 ) and also listed in 



Moustakas et al 



(|2010h : Col. (5) 



Physical size of adopte d aperture; Col. (6) Oxygen abundance information adopted from 



Moustakas et al. 



( 120101 ) : 1 = characteristic value for the galaxy is adopted; 2 = gradient is 
adopted and metallicity calculated for each aperture; 3 = L-Z relation derived value is adopted; 
Col. (7) The number of apertures selected from each galaxy; Col. (8) The sub-sample each 
galaxy belongs to, but the regions in the galaxy may belong to another sub-sample according 
to metallicity gradients, if available; L = low metallicity sample; M = intermediate metallicity 
sample; H = high metallicity sample. 
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Fig. 1. — Aperture and region selection of NGC5194 and NGC3627 on original Ha and 
24 /iin images (left and right panels), and NGC5194 on convolved Ha and 24 /im images 
(middle panel). The sizes are 7.05' x 11.25' (280x447 pc) for left and middle panels and 
3.83' x 6.15' (165x265 pc) for right panel. Circles are selected apertures and boxes are 

regions for local background subtraction. Central regions of both galaxies are not included 
as both are classified as having an AGN type nucleus. Red: 70 /im; green: 24 /im; blue: 
Ha. North is up and east is left. The two galaxies are chosen to show one example of a 

nearby extended galaxy with multiple resolved features and one example of a farther away, 
smaller galaxy with a smaller number of identified regions. 
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■4-1 

c 

D 
O 

o 




8.0 8.2 8.4 8.6 8.8 9.00.00 0.06 0.12 0.18 0.24 0.30 
12 + log(0/H) ct(z) 

Fig. 2. — The distribution of metallicity values for our data on the left panel and 
metallicity uncertainty on the right panel. The majority are in the intermediate metallicity 
sub-sample (8.4 < 12 + log(0/H) < 8.8) and have an uncertainty less than 0.15 dex. 
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38 39 40 41 

log[E(Ha) + 0.031E(24)] 



Fig. 3.— The 



( ICalzetti et al. 



SD o f 70 /im as a function of unbiased SFR indicator (surface density) 



20071 ) divided into three sub-samples according to the metallicity (Red 



squares, high metallicity sample with 12 + log(0/H) ^ 8.8; green pluses, intermediate 
metallicity sample with 8.4 < 12 + log(0/H) < 8.8; blue circles, low metallicity sample with 
12 + log(0/H) ^ 8.4). The typical error bar is shown as the plus on upper-left. The solid 
line is the best linear fit through the entire sample. The units are ergs • s _1 • kpc~ 2 for LSDs 

(luminosity surface density). 
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Fig. 4. — The residual between LSD of 70 /im and unbiased SFR indicator (surface 
density) as a function of the star formation rate density, SFRD (see caption of Figure |3] for 
color-code and symbol-code information). 1-a error is shown as the plus on upper- left for 
left panel. The units are ergs • s _1 • kpc -2 for LSDs and M ■ yr _1 ■ kpc~ 2 for SFRDs 

respectively. 
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CM 




8.0 8.2 8.4 8.6 8.8 9.0 - 1 .0 -0.5 0.0 

12 + log(0/H) log[L(24)/L(70)] 

Fig. 5. — The residual as a function of metallicity (left hand-side) and L(24)/L(70) ratio 
(right hand-side, see caption of Figure [3] for color-code and symbol-code information). The 
trend is clearly stronger for the residual as a function of metallicity, indicating that 
variations in the temperature of the warm dust (as traced by L(24)/L(70)) are not the 
dominant contributor to the scatter around the mean trend. The apparent correlation on 

the upper right part of the data on the right panel is due to the fact that the plot is 
actually 1/x versus x once IR emission dominates. Three short lines above the data points 
in the left panel show the mean metallicity uncertainties of each sample separately. The 

units are ergs • s _1 ■ kpc~ 2 for LSDs. 
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38 39 40 41 

log[S(Ha) + 0.031S(24)] 



Fig. 6. — Models of 3 different ages, 10 Myr, 100 Myr and 1 Gyr (blue dotted, green solid 
and red dashed respectively), with solar metallicity, overlaid on data (see caption of Figure 
[3] for color-code and symbol-code information). Typical 1-a error is shown as the plus on 
upper- left. Different ages of model do not show significant difference from each other. The 

units are ergs • s _1 • kpc~ 2 for LSDs. 



38 39 40 41 38 39 40 41 

log[S(Ha) + 0.031E(24)] 

Fig. 7. — Three models with different metallicities, 3.7 Z & (red dashed lines), 1.4Z Q (green 

solid line) and O.3Z (blue dotted line), with an age of 100 Myr, overlaid on data with 
corresponding metallicity sample (see caption of Figure [3] for color-code and symbol-code 
information) separately on each panel and the merged plot on the bottom right panel. The 
models reproduce the average trend of the data once the different metallicities are taken 
into account. Combined together, the data show the scatter which we attribute mainly to 
differences in metallicity. Typical 1-cr error is shown as the black plus on upper-left cornel. 

The units are ergs • s _1 ■ kpc~ 2 for LSDs. 
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Fig. 8. — The luminosity ratio L(24)/L(70) as a function of the star formation rate density, 

SFRD, for our sample regions (see caption of Figure [3] for color-code and symbol-code 
information). Three models with different metallicity (see caption of Figure [7] for line style 
information) are overlaid on data. The big filled symbols (blue circles for low metallicity 

sample, green diamonds for intermediate metallicity sample and red squares for high 
metallicity sample) represent the binned average (bins separated by vertical dashed lines) 
for each metallicity sample and the model lines follow reasonably well the average trend of 
the data. The units are M • s~ x • kpc~ 2 for SFRDs and ergs • s _1 • kpc -2 for LSDs separately. 
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Fig. 9. — The data at metallicity larger than 12 + log(0/H) = 8.4 are fitted with a line 
with fixed slope of unity in the £(70)-SFRD plane, 
log[E(70)] = log[SFRD] +43.028(±0.007) (68% envelop shown as dashed lines). The units 
are ergs • s _1 • kpc~ 2 for LSDs and M • yr -1 • kpc~ 2 for SFRDs separately. 
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Fig. 10. — On the upper panel, A(Ha) as a function of metallicity (see caption of Figure [3] 
for color-code and symbol-code information). The at tenuation of Ha is c alculated by using 



the L(Ha)+0.031L(24) as the intrinsic Ha emission (jCalzetti et al. 



20071 ). The attenuation 



increases with higher metallicity. The lower panel is the same plot as Figure [3J except the 

data is divided into three sub-samples depending on the A(Ha) values (red squares, 
A(Ha) > 1; green pluses, 0.25 < A (Ha) < 1; blue circles, A(Ha) < 0.25). However there 
may be objects with low extinction that have been excluded from the sample, due to our 
source selection criterion. Three short lines in the upper panel above the data points show 
the mean metallicity uncertainties of each sample separately. The units are 

ergs ■ s _1 ■ kpc -2 for LSDs. 



8.0 8.2 8.4 8.6 8.8 9.0 

12 + log(0/H) meon 



Fig. 11. — The calibration coefficient (scaled by a factor of 10~ 43 ) in Equation ([3]) derived 

from sub-samples with different (but overlapping) metallicity range as a function of the 
mean metallicity of the sub-sample. The error bar on mean metallicity shows the standard 
deviation for each sub-sample metallicity distribution and the error bar on the calibration 

coefficient shows the l-a dispersion of each sub-sample. The solid line is a linear fit 
through our six data points (circles). The filled circle is the mean calibration, adopted in 



the paper, for the intermediate and high metal 
calibration for integrated 70 /im from galaxies in 



icity samples. 



Calzetti et al 



he fi lled square is the 



fl2010h . The dashed lines 



indicate the calibration for Magellanic Cloud s (taking average of 12+log(0/H)=8.4 for 



Lawton et al. 



(120101 ) and the error bar is 



LMC and 12+log(O/H)=8.0 for SMC) from 
derived from their dispersion of 70 /zm/TIR ratio. These two points are not included in the 

fit shown in the Figure. 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 

log[L(24)/L(70) Whole] 



Fig. 12. — The L(24)/L(70) ratio summed up for all the regions in one galaxy as a 
function of the L(24)/L(70) ratio within the whole galaxy on the upper panel and the 
residual plot on the lower panel. Red circles are galaxies identified as having an AGN while 
blue squares are galaxies identified as not having an AGN. Open symbols identify galaxies 
with less than four regions each. For the AGN-free galaxies, the average dust temperature 
in HII regions is higher than the galaxy as a whole. 
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Fig. 13. — The residual between original and convolved (to the 70 /im PSF) photometry as 
a function of original Ha photometry, with local background subtraction. The convolved 
data are not corrected for any loss outside the aperture which has 16" radius. We use the 
high end 17 data points (~ 3%,blue points right to the dashed line) of the sample, i.e. the 
least contaminated apertures, to establish an initial constant aperture correction. The solid 
blue line indicates the mean of these 17 data points, which is also the adopted constant 
aperture correction value. The subscripts orig and conv here and in the paper are used to 
indicate the photometry from the original high resolution images and those from the 
convolved images respectively. The units are ergs ■ s -1 • kpc~ 2 for LSDs. 
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Fig. 14. — The residual between original and convolved photometry as a function of 
original Ha photometry, without local background subtraction for both. The low end 
shows an apparent flaring trend compared to the local background subtracted one, which 
suggests the need to apply the local background subtraction, especially for low luminosity 
data points. The units are ergs • s _1 • kpc~ 2 for LSDs. 
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Fig. 15. — The residual between original and convolved photometry as a function of the 

FWHM of the point source within the aperture. Filled circles are those with only one 
bright point source, contributing at least 50% of the flux within 2FWHM area, and open 
circles are those with one brighter point source but also a companion weaker point source, 
where the FWHM of the brighter source is measured. The dashed line is a linear fit to the 
filled circles. The aperture correction is not correlated with the compactness of the source, 
represented by the FWHM, within the apertures. The units are arcsec for FWHM and 

ergs • s _1 • kpc~ 2 for LSDs separately. 
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Fig. 16. — The residual as a function of constant aperture correction corrected photometry 
is fitted with a linear line to establish the surface brightness dependent aperture correction. 
The convolved photometry data have been corrected for the constant aperture correction 
and the subscript, c-corr, is used to indicate the constant aperture correction corrected 
photometry here and in the paper. The units are ergs • s _1 • kpc™ 2 for LSDs. 
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38.0 



38.5 



40.0 



40.5 



39.0 39.5 
log[S or;g (Ha)] 

Fig. 17. — The corrected photometry as a function of original Ha photometry on the 
upper panel, and the residual as a function of original Ha photometry on the lower panel, 

where the la dispersion (68% data envelope, dashed lines) is ~ 0.09 dex. The solid line 
indicates the one-to-one line. The subscript, 1-corr, is used to indicate the surface 

brightness dependent aperture correction corrected photometry here and in the paper, 
which is also the two-step aperture correction corrected photometry and thus the proxy to 
the correct photometry. In the rest of the paper, all luminosity and LSD data is from the 

1-corr photometry unless otherwise specified. The units are ergs • s _1 • kpc~ 2 for LSDs. 



38.0 38.5 39.0 39.5 40.0 40.5 

log[S orig (Ha)] 



Fig. 18. — The 1-corr photometry as a function of the original Ha photometry, while the 
aperture correction used here is established using only the brighter half, in LSD, of the 
data (circles right of the dashed line). The solid line indicates the one-to-one line and the 
dashed line separates the brighter half and the dimmer half. The squares show that the 
aperture correction established with high end half of data applies to the less luminous data 
as well, which suggests that our aperture corrections are robust. The units are 

ergs • s _1 • kpc~ 2 for LSDs. 



